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^! Abstract 

Oh: 

^ I The master formula approach to chiral symmetry breaking is used to analyze the phase 

shifts for TTTT scattering in the elastic region. The results are in excellent agreement with 
the data for the phase shifts up to the KK threshold. Our analysis shows that the vrvr 
data near threshold in the scalar channel favors a large quark condensate in the vacuum, 
i.e. (qq) ~ -(240 MeV)^. 



X 



SUNY-NTG-95-14 
May 5, 1995 



1. Introduction 

Chiral symmetry offers an important frameworlc for discussing hadronic processes at 
low energy ||l|. In the past, current algebra using the PCAC hypothesis [@] and chiral 
perturbation theory using a one loop expansion P] have been applied to describe processes 
involving pions. Both methods are limited to the threshold region because the former is 
an expansion around the soft pion point and the latter is an expansion around the chirally 
symmetric point. These points are unphysical in nature. Proposals to remediate these 
shortcomings based on dispersion techniques ||^, while well motivated from the point of 
view of unitarity, do not help in understanding quantitatively the role of chiral symmetry 
beyond threshold. 

Recently, two of us have proposed a general framework for discussing the role and 
consequences of chiral symmetry for processes involving pions at threshold and beyond 
1^. Our approach relies on a set of exact identities for on-shell pions, all of which following 
from a master equation for the extended S matrix. This equation, which has the structure 
of a reduction formula, captures the essence of chiral symmetry and generalizes the current 
algebra approach on-shell. Scattering amplitudes involving pions, nucleons and photons 
follow by reduction in terms of correlation functions, some of which are measurable. Our 
approach enforces chiral symmetry and unitarity to all orders in the pion momentum, and 
yields exact threshold theorems that are not necessarily contrived around the soft pion 
point or the chiral point. 

The purpose of this paper is to discuss in detail the tttt scattering amphtude as de- 
rived from the master formula in the elastic region and compare it to the available data. 
The amplitude is rewritten in terms of the scalar and vector pion form factors, as well as 
vacuum correlators. The occurrence of these terms is required by chiral symmetry and 
unitarity. The advantage of this approach is that not only can these quantities be esti- 
mated in the same fashion as chiral perturbation theory with fewer parameters, but they 
can also be approximated by well accepted forms suggested by experiment without fine 
tuning of the parameters or the theory. This can then be used to assess our understanding 
of the low energy physics without being hampered by the need to use a perturbation series 
or worry about unitarity. 

Given the importance of tttt scattering in hadronic physics, there are currently a variety 

1 



of theoretical approaches to vrvr scattering most of which allow for a good deal of phe- 
nomenological inspiration Since our purpose is to assess the constraints brought 

about by chiral symmetry alone, we will only compare our predictions with the data and 
one loop chiral perturbation theory |l^. For clarity, we will also display the one loop 
analysis following from the master equation by power counting in l/fj^ and minimality 
1^ . The latter involves only two parameters in the isospin symmetric limit whereas chiral 
perturbation theory requires four. Our results show that the master formula approach 
is able to reproduce the phase shifts quite well up to the KK threshold using educated 
inputs for the vacuum correlation functions. 

The formulas used throughout this paper are presented in the next section. Section 3 
applies these formulas to the one loop expansion in order to analyze tttt phase shifts. In 
section 4, the relevant form factors are modeled in the simplest form consistent with data 
and used with the master formula approach to fit phase shifts far above the threshold 
region. Comparison with chiral perturbation theory (xPT) results is made throughout. 
Our conclusions are summarized in section 5. 

2. The Calculation 

The master equations of the new formalism give a relation between the scattering 
amplitude and various correlation functions. Taking the case of tttt scattering, this rela- 
tionship is quite simple. Specifically 



using the Mandelstam variables s, t, and u. In terms of the pion vector form factor 
Fy, the pion scalar form factor F5, and the respective vector- vector and scalar-scalar 
correlation functions in the vacuum, Ily and XI^, the general result may be writtenQfor 
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The vacuum expectation value (a) was defined in to be 
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with m = 7 MeV being the average of the SU{2) quark masses. Taking a reasonable 
value for the quark condensate (qq) = — (240MeV)'^ gives (a) ~ —40 MeV. We recall that 



(a) = corresponds to the Gell-Mann-Oakes-Renner (GOR) relation [|r^. A discussion 
of the overall dependence of our result on (a) is given at the end of section 4. 

The tree amplitude was first calculated by Weinberg [||] and must appear in any low 
energy theory. Also appearing is the vector and scalar expressions showing in a model 



independent way the importance of the p contribution |]T3| and the a state consisting of 
the two scattered pions [|T^. The contribution from the "rest" is a correlation function 
of four one-pion reduced axial vectors . Since this is not a measured quantity, its one 
loop form is quoted. Specifically, 
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with 

-^'^ > = - f - ^arccoth^l - (3) 

representing the finite part of the pion-propagated loop. 

Equation (|]) is an exact result. However, there is no known way at present to funda- 
mentally calculate the correlation functions or form factors of this result. Instead, either a 
loop expansion or experimental fit must be used. For both methods, the different isospin 
contributions can be extracted from the pion scattering amplitude in the usual manner 



T°(s,t) = 3A{s,t,u) + A{t,u,s) + A{u,s,t) 
T^(s,t) = A{t,u,s) - A{u,s,t) 
T^is,t) = A{t,u,s) + A{u,s,t) 



T\s,t) = 3271 J y(2l+l)Pi(cose)e'^'^'hm6Us). 

\ s-Aml f^^ 

Here 6 is scattering angle in the center of mass frame. Since elastic scattering will always 
be considered, no part of the amplitude will be lost to another channel and so the phase 
shifts will always be real. 

To obtain the phase shift, 6{ , an integration of the respective over the Legendre 
polynomial with index I must be performed. This leads to a unitarity condition on the 
scattering amphtude 



fi is) = yi — ~ ~6A^ {s,cose)Pi{cose) 

Im// = |//f (4) 

which serves as a consistency check on the perturbative expansion employed. 
Two of the various ways to express the phase shift are 

Sf{s) = arccot (^^) = ^ arcsin (2 Re //(.)) (5) 

The first expression is the definition for the phase shift and is therefore preferred. However, 
it is sometimes better to use the second expression as will be expounded upon below. The 
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equivalence of these two expressions for the phase shift may also serve as a measure of 
unitarity. 



3. Loop expansion 

Calculation of the amplitude in loops is really an expansion in powers of s/(47r/^)^. 
This expansion has relevant corrections only for ^/s > 0.5 GeV. The amplitude is purely 
real at tree level and only at one loop is there any contribution to the imaginary part. It 
is therefore easy to see that such an expansion does not preserve unitarity (eq. (^). For 
the same reason, the arcsin rather than the arccotangent is a better approximation to the 
true phase shift at one loop since two terms of the expansion are known for the real part 
as opposed to only one term for the imaginary part. Therefore the arcsin will be used for 
the graphs of this section. 

It is not important to quote the one loop form of the scattering amplitude since it is 
worked out in |^ except to note that the function J{q^) given by equation @ appears in all 
the form factors. In the t and u channels, the argument of J will have a cos 9 dependence. 
The integral over the Legendre polynomials then must be integrated numerically. 

As explained in 0, only two constants appear in the amplitude to one loop: ci and ci. 
There it was shown that Ci = ny(0) = 2fl'F'y{Q) = 0.031 ± 0.001 to one loop using the 
measured value for the vector charge radius {r'^)v = 6Fy(0) = 0.42±0.02 fm^ [Q. In f^, 
various scattering lengths were calculated and compared with experiment to give seven 
determinations of ci. The weighted average of these values is IGvr^ci = 2.46 ±0.55. Table 
1 shows how the scattering lengths and range parameters in the master formula approach 
using the middle value of ci compare with experiment, the tree calculation, and xPT0 
(using arcsin for all theoretical values). The phase shifts 6q, 6l, and 6q are plotted as the 
solid line in figures 1,2, and 3. All curves in these figures cannot be extended beyond 45 
degrees because the one loop result keeps growing like In pushing the argument of 
the arcsin above one. 

In all three graphs, both results fit the data well up to about 450 MeV. The subsequent 
deviation is most likely an artifact of the loop expansion, but a two loop calculation is 



needed to clear up the ambiguities [T^. The shaded region for the master equations in 



^We used the central values of li = —0.6, I2 = 6.3, Z3 = 2.9, and I4 = 4.3 as quoted in 



6q shows the variation due to the uncertainty in ci. This shows that up to the point of 
deviation, the lack of accuracy of Ci does not matter much. Therefore the middle value 
for Ci was used for all other calculations in this paper. The master formula approach 
is closer to the experimental value for 6l and 6^ even though it has half the number of 
tunable parameters compared to xPT. For 6q the master formula result is too low at 
threshold, but the overall shape leads more towards the higher energy values before the 
dependence dominates. Since the one loop results in the master formula follow from 
a 1//^ expansion and minimality the present results show the overall consistency of 



the approach. Minimality is at the origin of KSFR |T2 . 

Taking the arccotangent form for the phase shift to one loop allows the result to be 
in the full range of to 180 degrees. It was argued above that the imaginary part will 
be too small since only the first term of the perturbation series is present to one loop. 
This would result in a lower phase shift than the full result of the theory should give. 
Indeed this turns out to be the case as seen by the dashed line in figures 4, 5, and 6. This 
is especially apparent for 6l which shows no sign of a p resonance from this analysis for 
either the master formula or xPT since it is impossible to see resonances from a strict 
perturbative expansion due to violation of unitarity. 

4. p and a saturation 

The one loop result employed above works well near threshold but, as was seen, cannot 
model even the simplest resonances. For this reason an alternative approach would be 
preferred. The master formula gives a general result which is not dependent on the loop 
expansion and so allows such an alternative. Modeling the form factors and correlation 
functions that appear in eq. with reasonable forms suggested by experiment should 
lead to an accurate result without much fine tuning of the parameters or of the theory. 

The only problem is that direct experimental measurements do not exist for some of 
the correlation functions needed. At least for the vector channel it is well accepted that 



the features are dominated by the p and a good parameterization is |T^, |Tj] 



^v{q') = — / .\ . (6) 



with mp = 770 MeV and a momentum dependent width 
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For vector dominance 7 is zero. Fitting the curve to the data for Fy requires 7 = 0.3 and 
is plotted in figure 7. This corresponds to {r'^)v = 0.51 fm^ which is reasonably close to 
the experimental value quoted above. 

The vector- vector correlation function, ny(g^), can be assumed to be dominated by 
the p resonance as well. Therefore the imaginary part may be modeled by a delta function 
for the p and a theta function so that the perturbative QCD result is recovered at large 
energies 

lmUy{s) = fl5{s - ml) + ^e{s - s,) 
^ ' 47r 



with So — 1.2 GeV as a reasonable value taken for the continuum threshold |T8[. The 
real part may be obtained by a once subtracted dispersion relation. The same can be 
done for the scalar-scalar correlation function 115(5^) with a twice subtracted dispersion 
relation due to the form of the perturbative QCD result. Smearing out the delta function 
by giving a width to the resonances, yields 



n.w = ^-iL_^_J_in(i_±)+nr*-(0) 

mj^ — s — impTp 4:71^ \ sqJ 

L^s{s) = — 2 ■■ — A 2 4 ^ 1^ U 

m^ — s — tmiL 1 47r^m^ V Sq 



+fjur'-{o)+sur'- '(0) 



Additional resonances can be added to the scalar channel if needed by adding further terms 
similar to the first term in the equation for II^. Vacuum dominance in the vector channel 
requires ny^^^'(O) = 0, however this requirement can be relaxed for scalar particles and 
used to enable a better fit to data. Therefore let a = /^n5^^^'(0) and /? = f^H^^^^' '(0) 
which are both dimensionless. The widths are zero for spacelike momenta and are taken 
to be eq. for the vector width and 

for the scalar. 
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Expanding the resonance forms in s and matching coefficients with those of the one 
loop result in the vector channel gives fp = (140 ± 4) MeV in agreement with the 
experimental value of 144 MeV. Furthermore, if one resonance is assumed in the scalar 
channel, matching threshold parameters for the scalar channel gives a = —2 x 10~^, 
(3 = —5 X 10~^, and fi = 125 MeV. These values were used for the fits to the phase shifts 
and are close to those used in [0 (which in our notation would be: mi = 540 MeV, 
/i = 90 MeV, and So = 1.1 GeV). 

Finally considering the scalar form factor, Fs{q^), there is not much to hint at its form. 
At tree level F5(0) = — I/Ztt, so a reasonable form assuming the same one resonance 
as used for lis is 

1 hiinl 
^sis) = - — —2 



/tt nil — s — irriiTi 



Using the empirical determination for the scalar charge radius 6Fs(0)/F5'(0) = 0.7±0.2 
fm^ 1^ and taking hi = 1 gives mi = 580 MeV. Fitting to the data for 6q requires the 
only free parameter, Fi, to be 175 MeV. Figure 7 shows that |Fy| and Ily both compare 
well with the data |1^. The imaginary part of F5 normalized by F5'(0) is comparable to 
a dispersive analysis shown by the dashed line. 

Even though there is no definite resonance below the /o(980), many have interpreted 
the large enhancement of the 6q data in the region of ^/s = 500 — 800 MeV (see figure 
4) to be due to a cr "particle" with mass around 500 MeV and width ^ 700 MeV |T3 



This scenario is consistent with our value for mi but our Fi is much smaller. It should 
be noted that in the master formula approach a large contribution from the t-channel p 
resonance is also seen. This is in agreement with conjectures set forth by other authors 
although it is not sufficient to alone fit the data in this region. Figure 4 shows that 
(mi,ri) = (580, 175) MeV reproduces the data quite well up to about 650 MeV. Above 
this the ffi rises slower than the data. 

Enhancement of the 6q resonance at higher energies and also a leveling of the sharp 
peak in Us (figure 7) can be done through the addition of further resonances. It was 
found that including both of the next two 0^(0++) resonances: /o(980) and /o(1300) 
worked well. The ffi to the data is shown by the solid line in figure 8 and the form factors 
in figure 9. Values for the constants employed are shown in table 2. The result strongly 
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depends on the values of the /ij's and only weakly depends on the /j's. Since the first 
resonance at 560 MeV fit the data very well, its constants were retained. For such a 
simple model of the resonances the fit is quite good up to the KK threshold. Here the 
strange particles' contribution and the role of the "rest" term, two problems which may 
be related, need to be taken into account and will be studied elsewhere . This fit shows 
in Yls (see figure 9) that even a narrow resonance around 580 MeV can be obscured by 
wide neighboring resonances. Further analysis is required. 

The results are very good in the vector channel as well, giving the p resonance as 
expected in figure 5. This is primarily due to three reasons: there is a large amount of data 
available in this channel and so the model used for the form factors is phenomenologically 
based, there is only one resonance of importance in the energy region of interest, and the 
"rest" does not contribute much to one loop showing this is probably the case in general. 
In addition the scalar part of the amplitude has only a small contribution and hence the 
one resonance or the three resonance scalar fit discussed above give identical results. For 
definiteness, one scalar resonance was used. 

The "rest" term is important to the result. In fact it is the only term which 
significantly contributes to the imaginary part of the scattering amplitude. This term 
was not taken into account in the resonance forms and so the phase shift is nearly always 
0° as seen by the solid line in figure 6. However, using the one loop expression for the "rest" 
is justified up to about 0.5 GeV as seen in the last section. Using the one resonance fit and 
including the one loop "rest" agrees very well with data and is also plotted in figure 6 as 
the dashed line. The imaginary part of the expanded "rest" dominates at high momenta 
due to the dependence and hence the result begins to deviate around 450 MeV. Also 
figure 6 includes the arccotangent form of the xPT result given by the dotted-dashed line. 

To see the threshold behavior, the combination — is plotted for the master formula 
approach using one resonance and for both the master formula approach and xPT in the 
one loop approximation in figure 10. All three give comparable results near threshold 
where there is data. This combination is no longer as important a parameter to plot since 
detailed information on both phase shifts independently is now available. 

As discussed above, the value of (a) = —40 MeV was used in all calculations. However, 
the GOR relation gives (a) = 0. Furthermore, if the quark condensate (gg) is zero then 
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(a) = —93 MeV. The dependence of the master formula solutions on the value of (a) can 
be seen best in the scalar channel which contains the largest contribution. Only the real 
part of the scattering amplitude is affected. Figure 11 shows the dependence of 6q on 
these three values of (a) for the one loop master formula calculation. This graph shows 
that (a) = —93 MeV, corresponding to a zero quark condensate, has the wrong shape at 
threshold. The other two values for (a) are comparable with the data. Figure 12 shows 
how the combination 6q — b\ changes with the three different values of (a). Here the 
—40 MeV value seems to fit best. 

As was pointed out by [^, the difference — mass (498 MeV) is a direct 



measurement of CP violation. The most recent data [21] gives 29.2 ±3° from np ttttA. 



A more direct measurement from K ^ tttt decay is not available, although analyses based 



on chiral perturbation theory have been carried out A weighted average of the data 
prior to 1975 gives ^3] 41.4° ±8.1°. Using the master formula resonance result (with one 
scalar resonance) gives 29.8°. Use of the one resonance result with the one loop "rest" 
added is not reliable at this energy but gives 53.0°. The one loop expansion in general can 
only be used in the arccotangent form since the arcsin form of Sq does not extend to the 
mass. This gives 37.2°. In comparison, the arccotangent form for xPT gives 48.3° (to 
be specific, quote 45° ± 6° where they use the arccotangent with the real part only). 

5. Conclusions 



The master formula approach to tttt scattering offers powerful constraints on the scat- 
tering amplitude solely on the basis of chiral symmetry and unitarity. As a loop expansion, 
it fits the threshold region as well as chiral perturbation theory with half as many param- 
eters to fix, and is compatible with KSFR. As a resonance fitting of the form factors and 
the vacuum correlation functions, it is in excellent agreement with the data up to the KK 
threshold. We stress that the various form factors and correlation functions are amenable 
to fundamental estimates in QCD or empirical measurements. With a minimum number 
of assumptions, this approach provides a testing ground for an empirical understanding 
of scalar and vector correlators from the vrvr data. Also, the tttt data near threshold in 
the scalar channel, although not very accurate, suggests that the vacuum supports a large 
quark condensate, and points in favor of the GOR relation. Better measurements at 
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threshold will settle these important theoretical issues. 

The present results also show that chiral symmetry constraints go far beyond the 
threshold region. In particular, we have found that in the vector channel the p is domi- 
nant, while in the scalar channel the presence of a low-lying though not necessarily broad 
resonance is unavoidable. Since our approach shows clearly the interplay between cor- 
relations in the vacuum and measured scattering amplitudes, it is well suited for lattice 
estimations. Improvements can even be made to the master formula approach by an ex- 
pansion to S't/(3) and inclusion of kinematical nucleons. These additions should allow a 
much better understanding of the inelastic region and will be carried out elsewhere |]19|. 
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igure Captions 



Fig.l. The 6q phase shift using the arcsin for the one loop expansion. The shaded 
region is the master formula result taking into account the uncertainty in ci. The 
central value of IGvr^ci = 2.46 is given by the middle solid line. The tree result is 
shown by the dashed line and the dotted-dashed line is the xPT result. Data is 



from |24| below 400 MeV and 25 otherwise 



Fig. 2. The 6l phase shift using the arcsin for the one loop expansion. The solid line 
is the master formula results, dotted line is the tree result, and the dotted-dashed 
line is the xPT result. Data is from p4, 2^ . 



Fig. 3. The Sq phase shift using the arcsin for the master formula (solid line), tree 



(dashed line), and xPT (dotted-dashed line). Data is from p8 



Fig. 4. 5q using the arccotangent. Here the solid line is the one resonance fit of the 
master equations for one resonance, the dotted-dashed line is xPT to one loop, and 
the dashed line is the master formula result to one loop. Data is from p^, |25| . 

Fig. 5. 6l using the arccotangent for the master formula with one resonance (solid 
line), master formula expanded to one loop (dashed line), and xPT (dotted-dashed 
line). Data is from |23, 25 . 



Fig. 6. Sq using the arccotangent for the master formula approach using one reso- 
nance (solid line) and the one loop expanded "rest" along with the one resonance 
(dotted line). The xPT result using arccotangent (dashed line) is also shown. Data 
is from p8i. 



Fig. 7. The form factors for the one resonance fit of the master formula approach. 
Data for Fy and Ily is from [|T^ (see for details) and the estimate for from 
a dispersive analysis 0. 

Fig. 8. Fitting the 6q to three resonances in the master equation approach. Also 
plotted are the master formula to one loop (dashed line) and xPT (dotted-dashed 



line). Data is from ||2J, ^ 
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Fig. 9. The form factors for the three resonance fit to the master formula (sohd hne). 
Other data is as in figure 7. 

Fig. 10. Sq — Si with the one resonance fit for the master formula result (solid line), 
the one loop master formula result (dashed line), and xPT (dotted-dashed hne). 
Data is from p4| . 

Fig. 11. The 5g phase shift for the one loop expansion of the master formula using 
{a) = (solid line), -40 MeV (dotted line), and -93 MeV (dashed line). These 
values correspond to {qq) = — (290MeV)^, —(240 Me V)^, and respectively. Data 
is from Q ^. 

Fig. 12. 6q —6l for the master formula one resonance fit using (a) = —93 MeV (solid 
line), —40 MeV (dashed line), and (dotted-dashed line). These values correspond 
to {qq) = 0, -(240MeV)^ and -(290 MeV)^ respectively. Data is from [||]. 
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Table 1: Comparison of the experimental scattering lengths and range parameters to those 
predicted by the tree, xPT, and one loop master formula calculations. For the master 
formula, the values ci = 0.031 and IGvr^ci = 2.46 were used. Data is from |^ and xPT 
results are from ^. 





Experiment 


Tree 




Master 


"o 


0.26 ±0.05 


0.16 


0.20 


0.20 


ml Wo 


0.25 ±0.03 


0.18 


0.26 


0.24 


2 1 


0.038 ±0.002 


0.030 


0.033 


0.038 


al 


-0.019 ±0.021 


-0.045 


-0.041 


-0.041 


ml bl 


-0.082 ±0.008 


-0.089 


-0.065 


-0.068 



Table 2: The parameters used for the three resonance fit of the master formal approach. 
All constants are in MeV except for a, /3, 7, and the /ij's which are dimensionless. 



mi 


580 


/i 


125 




175 


m2 


980 


/2 


300 


r2 


400 


ms 


1300 


/3 


300 




400 


m 


7 


hi 


1 


(^) 


-40 




1200 




13 


a, 7 


-0.17, 0.3 


fp 


140 


h 


-11 




-1.6 X 10-3 
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